Intrinsic vs. extrinsic anomalous Hall effect in ferromagnets 
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A unified theory of the anomalous Hall effect (AHE) is presented for multi-band ferromagnetic 
metallic systems with dilute impurities. In the clean limit, the AHE is mostly due to the extrinsic 
skew-scattering. When the Fermi level is located around anti-crossing of band dispersions split by 
spin-orbit interaction, the intrinsic AHE to be calculated ab initio is resonantly enhanced by its non- 
perturbative nature, revealing the extrinsic-to-intrinsic crossover which occurs when the relaxation 
rate is comparable to the spin-orbit interaction energy. 
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Early experimental works on the Hall effect in ferro- 
magnetic metals led a semi-empirical relation of the Hall 
resistivity p xy to a weak applied magnetic field H and 
the spontaneous magnetization M both along the z di- 
rection; p xy = RhH + AttR s M with the normal and the 
anomalous Hall coefficients Rh and R s , respectively 0. 
This anomalous Hall effect (AHE) has been one of 
the most fundamental and intriguing but controversial 
issues in condensed-matter physics 0, IE H IE IE • ^ 
has not been clarified yet if the AHE is originated purely 
from extrinsic scattering or has an intrinsic contribu- 
tion from the electronic band structure, which penetrates 
even recent debates on the interpretation of the experi- 
ments |E |E El • Theoretically, a unified description of 
both intrinsic and extrinsic contributions is called for but 
has not been considered seriously. It even reveals their 
nontrivial interplay and crossover and explains the AHE 
in a whole region from the clean limit to the hopping- 
conduction region (see Fig. 0J, which are main goals of 
the present study. 

Dissipationless and topological nature of the Hall effect 
has been highlighted by the discovery of the quantum 
Hall effect in two-dimensional electron systems under 
a strong magnetic field. In the Stfeda formula 0] of the 
electric conductivity tensor = a 1 ^ + a{- , in ideal 
cases where the Fermi level is located within the energy 
gap, the Fermi-surface contribution <7y vanishes and the 
quantum contribution <r/ 7 yields the TKNN formula |l3j 
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with the electronic charge e, the Planck constant h = 
2irh, the Fermi distribution function /(e), and the anti- 
symmetric tensor e^. We have introduced the eigenen- 
ergy e n (p), the Berry curvature b n (p) = V p x a n {p) and 
the Berry connection a n {p) = i(n,p\ V p |n,p) of the gen- 
eralized Bloch wave function \n, p) with the band index n 



and the Bloch momentum p. Each band has a topological 
integer called the Chern number C n = — J n^n^ nip) ■ 
Their summation over the occupied bands determines the 
integer v (Chern number) for the quantization ct*°* = 
ve 2 /h in insulators. Then, adiabatic semi-classical wave- 
packet equations have been devised to incorporate this 
Berry-curvature effect into the equations of motion |14j . 

Historically, Karplus-Luttinger initiated an intrinsic 
mechanism of the AHE in a band model for ferromagnetic 
metals with the spin-orbit interaction, which coincides 
with the TKNN formula a TI<NN HE 111- This reflects 
that the spin-orbit interaction bears a nontrivial topolog- 
ical structure in the Bloch wave functions of ferromagnets 
by splitting band dispersions, which originally cross at a 
certain momentum po, with a transfer of Chern num- 
bers among the bands. This phenomenon called "parity 
anomaly" has a non-perturbative nature |l7j ] , and points 
to an importance of the anti-crossing points with a small 
gap 2Ao, which is identified with the spin-orbit interac- 
tion energy eso- When the Fermi level is located around 
such anti-crossing of dispersions, as found in recent ab 
inito calculations for SrRu03 an( l t ne bec Fe [l9| . 
the magnitude of o~ T ? NN is resonantly enhanced and ap- 



proaches e /h 
e 2 /ha - 10 3 n 



- 3.87 x 10 in two dimensions and 

1 cm -1 in three dimensions with the lat- 
tice constant a w 4 A 0, 0] . 

On the other hand, adiabatic semi-classical Boltzmann 
transport analyses suggest that impurity scattering pro- 
duces the AHE through the "skewness" [E IE or the 
side jump 0, 0. The skew-scattering contribution di- 
verges in the clean limit as cr|y CW = a xx S = j^^jf-S 
with the life time r and the Fermi energy Ep. Here, 
S ~ ssoVi m p/W 2 (<^i 1) is the skewness factor with W be- 
ing the bandwidth or the inverse of the density of states 
and Wimp the impurity potential strength. 

Generic model that fully takes into account both the 
"parity anomaly" associated with the anti-crossing of 
band dispersions and the impurity scattering can be ob- 
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FIG. 1: (Color online) Energy dispersions for the Hamiltonian 
H in Eq. (0. 



tained by expanding the Hamiltonian at a fixed p z with 
respect to the momentum p measured from the originally 
crossing point po of two dispersions; 



H +H [mp = -A a z +vp-d-xe z 



P 

2m 



T-imp 

with the position r of electron, the Pauli matrices & = 
(a x , a y ,a z ), the unit vector e z in the z direction, and an 
impurity at a position r; mp . The first term corresponds to 
the level splitting 2Ao = £so of two bands at the anti- 
crossing momentum. The second term gives the linear 
dispersion with the velocity v. The third term represents 
the quadratic dispersion with an effective mass m, whose 
anisotropy has been neglected since it is unimportant. 
This model has two band dispersions shown in 

Fig. n] Henceforth, the bottom of the lower band is cho- 
sen as the origin of the energy and the bottom of the 
upper band denoted as E ies = £-(0) is taken as an en- 
ergy unit. The model possesses the gauge flux b± — 
TA 2 A /2A3 with A p = i/AV + A 2 [U EH. In the 
case of resonance, Ef £ [E les — 2Ao, -Ei-es], ^ y NN ap- 
proaches the maximum value e 2 /2ft. Away from this res- 
onance, dominant contributions from the momentum re- 
gion around p = cancel out each other or do not appear, 
leading to a suppression of a^ NN '(« (e 2 /h)(sso/ Ef)), 
where the perturbation expansion in Sso is justified. 
Therefore, the present model, Eq. ©, can be regarded as 
a generic continuum model for a momentum region that 
gives a major contribution to the AHE. By definition of 
the anti-crossing, A does not change its sign as a func- 
tion of p z - This removes a concern that the integration 
over p z might lead to a cancellation. 

We employ the Keldysh formalism for non-equilibrium 
Green's functions, which has recently been reformulated 
for generic multi-component systems We consider 

Green's functions and self-energies under an applied con- 
stant electric field E = (0,E y ); G a (s,p) and E"(e) with 
a = R, A, < for the retarded, advanced and lesser com- 
ponents, respectively, e and p represent the covariant en- 
ergy and momentum in the Wigner representation |22| . 



G a (e,p) and S a (e) can be expanded in E y as 

G a {e,p) = G%(e,p) + ehE y G a Ey (e,p) + 0(E 2 y ), (3) 



X<*(e) = S«(e) + eHE y E%(e) + 0(£ 2 ) 



(4) 



Henceforth, functionals with the subscripts and E y de- 
note those in the absence of and the gauge-invariant lin- 
ear response to E y . Due to the 6- functional form of the 
impurity potential, the self-energies are local. Gq ' A sat- 
isfies the familiar Dyson equation, 



G*' A (e,p) = [e-Ho(p)-£Z' A (e)}-\ 



(5) 



R A <L 

The self-consistent equations for G E ' ' are obtained by 



expanding the Dyson equation in the electric field [2/ 



It is convenient to decompose G E and into two; 



G< y (e,p) = G< y>I (e,p)d e f(e) + G< ytII (e,p)f(e$) 

£< y (e) = t< yJ (e)d e f(e) + t< yJI (e)f(e) (7) 

Gi ytII {e,p) = G A y (e,p)-G%(s,p) (8) 

= t%{e)-±%{e). (9) 

G E j and / can be self-consistently determined from 
the quantum Boltzmann equation in the first order in E y , 



'EyJ^O 



— ^Eyj^a n E v j 
+ l - {{t A - t R ){d Py G A ) + {d Py G R )(t A E*))l0) 



fiy> G A - G R 



with the velocity Vi(p) = d Pi H (p), while G R ' A and T, E 
are determined from the other self-consistent equation 



R,A 



AR,A _ A,R,A^R,A/^R.A 

Ey - ^0 Ey U Q 



G*> A v v (d e G*' A ) - (d £ G R ' A )v y G^ A ) . (11) 



R A 

We can exactly calculate the self-energies E ' and 
up to the riimp-linear terms by means of the T- 



Sfl,A,< 
J E u 



matrix approximation; 

= n imp f*' A (e) 



f R '%) 



(12) 



for the zeroth-order in E y and 



d 2 p 



£ Ey ,M = ^ P f R (e)l—^G< y>I (e,p)T A (e) (14) 

:G R E A (e,p)trm) 



^■eI A ( £ ) = n imp f R ' A (e) 



for the first-order in E y . We solve Eqs. © , ifT^ and (|T3|) 

self-consistently to obtain Gq' A and T, R ' A . Next, they 



FIG. 2: (Color) (a) The total anomalous Hall conductivity 
a^y against Ef and h/r in an energy unit of E TeB — 1.0. (b) 
The intrinsic contribution cr^y for the same parameter values. 
Note the difference of the scales for a xy in (a) and (b). 
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are plugged into Eqs. (frU|) and (fH|) to solve G^ j and 



£< r self-consistently. G* A and ±f A are obtained from FIG - 3; ( Color online ) CT 4 and °*v as a function of ft/r for 



Eqs. {TTJ and {TSJ| , and hence G^ 7J through Eq. |(8J|. 
The conductivity tensors are calculated from 



the same parameter values as Fig. with Ef 
1.5 for (a), (b), and (c), respectively. 
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0i(p)Gl,„(e,p) MP) 



with i, j = x, y. Note that Eqs. (|16|l and l|17(l are along 
the same spirit as the Stfeda formula 12]: this approach 
provides the diagrammatic treatment for the Stfeda for- 
mula 22]. Effects of the vertex corrections to Gg. u 

cancel each other, and hence we can regard ajy as an 
intrinsic contribution. By contrast, the Fermi-surface 
contribution erL, suffers from a vertex correction. While 
the intra-band matrix elements correspond to the con- 
ventional description of both a xx and a xy based on the 
scattering events, the inter-band ones contain a finite in- 
trinsic contribution to the AHE as a part of the Berry- 
curvature term and is generically expressed as 



/ int 



dp 



{2nh) 



(e„(p) -£ n '(p)) 



xdJ(e n {p))Im((np\V p \n'p) x {n'p\ V p \np)) e .(18) 

Figure Et a ) shows the numerical results on a x °y — 
<7 xy + <r X y against the Fermi energy Ef and the Born 
scattering amplitude K/t = ni mp vf mp m for a typical set 
of parameters, A = 0.1, 2mv; mp = 0.6, 2m\ 2 = 3.59, 



and the energy cutoff is taken as E c = 3.0 in an energy 
unit of E rcs = 1.0 H3- In the clean limit H/t -> 0, a x °* 
diverges in accordance with the skew-scattering scenario. 
Strength of the divergence is proportional to Ef in the 
low electron-density limit, and the sign is inverted around 
Ef = e+(0) = E xes — 2Aq. Sign of the skew-scattering 
contribution also changes by the sign change of v- lmp . Fig- 
ure Etb) shows the intrinsic contribution ct™* calculated 



by imposing 2 



R,A,E 
E*. 



= for the same set of parameters. 



Under the resonant condition for Ef, a™* becomes of 
the order of e 2 /2h. With increasing h/r, it gradually 
decreases solely due to damping of quasi-particles. 

By contrast, with increasing h/r, the extrinsic skew- 
scattering contribution rapidly decays (Fig. reflect- 
ing that it originates purely from intra-band processes 
and hence the skewness factor S remains of the order of 
ssoVi mp /W 2 . In moderately dirty cases, the total con- 
ductivity nearly merges into the intrinsic value. Namely, 
there appears a crossover from the extrinsic regime to 
the intrinsic as h/r increases. Especially in the reso- 
nant case (Ef — E Tes ) shown in Fig. 01(b), the intrinsic 
contribution is significantly enhanced and the extrinsic- 
to-intrinsic crossover occurs at H/t ~ eso- For a small 
ratio of eso/Ef ~ 10~ 3 -10~ 2 dominance of the 

intirnsic AHE is realized within the usual clean metal. In 
reality, the total Hall conductivity is the sum of the con- 
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FIG. 4: (Color online) Scaling plot of a xy versus a xx for the 
same sets of parameter values as in Fig.^Jb) except 2mvi mp . 



tributions from all over the Brillouin zone. Since skew- 
scattering contributions from other momentum regions 
are always subject to a similar rapid decay, the above 
extrinsic-to-intrinsic crossover still occurs unless contri- 
butions from all the anti-crossing regions of band disper- 
sions are mutually canceled out by accident. 

Figure 01 shows the logarithmic plot of a xy against a xx 
for the same set of parameters as Fig.^Jb) except for the 
impurity scattering strength V{ mp , which is varied for dif- 
ferent curves. In the clean limit, the curves nicely follow 

&xy ^ @ xx 

and the ratio a xy /a xx is proportional to «i mp 
for a fixed r. As a xx = 2{e 2 /h)Epr /h decreases with de- 
crease in r, the relation exhibits a upward deviation from 
the linear one, signalling the crossover to the intrinsic 
regime. A smaller impurity potential strength Vi mp en- 
larges the region of the constant behaviour of a xy . (Note 
that we change ni mp to control h/r.) Careful experiments 
are required to test the prediction of the crossover at 
low temperatures. The magnitude of a xy in the intrinsic 
regime is consistent with experimentally observed values 
a xy = 10 2 — 10 3 J7~ 1 cm _1 in a cr^y-constant region of Fe- 
and Ni-based dilute alloys 0, and SrRu03 and metallic 
foils [25j . Further decrease of r again changes the scal- 
ing behaviour to a xy oc a xx , which almost agrees with 
recent experiments on a disordered pyrochlore ferromag- 
net Nd 2 (Moi_ x Nb x ) 2 07 26] and on Lai-^Sr^CoOa |25j . 
This exponent approximates to the value expected for the 
normal Hall effect in the hopping conduction regime |27| . 

Now the source of the confusion over decades is clear. 
The skew-scattering contribution, though it is rather 
sensitive to details of the impurity potential and band 
structure, can be larger than e 2 /h in the superclean 
case h/r <C eso, but decays for eso < h/ T - The 
side-jump contribution is smaller and of the order of 
(e 2 /h)(esol 'Ep) 7\. Therefore, the intrinsic one, which 
is of the order of e 2 /h under the resonant condition, is 
dominant over a wide range of the scattering strength 



eso < h/r < Ep (clean case). Although Luttinger recon- 
sidered the Karplus-Luttinger theory |2| and gave an ex- 
pansion of a xy in Wimp, including the skew-scattering con- 
tribution as well 5] , it fails to reveal the above crossover 
in the space of Ep, eso and fi/V. 

In conclusions, we have shown that the AHE is 
determined by the intrinsic mechanism when (i) the 
Fermi level is located around an anti-crossing of band 
dispersions in the momentum space, (ii) consequently 
the magnitude of a xy is comparable to e 2 /(ha) ~ 
10 3 f2 _1 cm -1 , and (iii) the resistivity p xx is larger than 
(ha/e 2 )(eso/Ep) ~ 1-10 fitl cm. With these conditions, 
first-principle calculation can give an accurate prediction 
of a xy . The present work resolves a long standing con- 
troversy on the mechanism of the AHE in a whole region. 
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